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ABSTRACT 

We reconsider the enhancement of primordial gravitational waves that 
arises from a quantum gravitational model of inflation. A distinctive feature 
of this model is that the end of inflation witnesses a brief phase during which 
the Hubble parameter oscillates in sign, changing the usual Hubble friction 
to anti-friction. An earlier analysis of this model was based on numerically 
evolving the graviton mode functions after guessing their initial conditions 
near the end of inflation. The current study is based on an equation which 
directly evolves the normalized square of the magnitude. We are also able 
to make a very reliable estimate for the initial condition using a rapidly 
converging expansion for the sub-horizon regime. Results are obtained for 
the energy density per logarithmic wave number as a fraction of the critical 
density. These results exhibit how the enhanced signal depends upon the 
number of oscillatory periods; they also show the resonant effects associated 
with particular wave numbers. 
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1 Introduction 



The case for a phase of accelerated expansion (inflation) during the very 
early universe is strong. One reason is that we can observe widely separated 
parts of the early universe which seem to be in thermal equilibrium with 
one another [I]. If one assumes the universe never underwent a period of 
inflation, there would not have been time for this thermal equilibrium to be 
established by causal processes. Without primordial inflation the number of 
causally distinct regions in our past light-cone at the time of recombination 
is over 10 3 , and it would be 10 9 at the time of nucleosynthesis. 

There is no strong indication for what caused primordial inflation. A 
natural mechanism for inflation can be found within gravitation - which, 
after all, plays the dominant role in shaping cosmological evolution - by 
supposing that the bare cosmological constant A is not unnaturally small 
but rather large and positive. Because A is constant in space, no special 
initial condition is needed to start inflation. We also dispense with the need 
to employ a new, otherwise undetected scalar field. However, A is constant 
in time as well, and classical physics can offer no natural mechanism for 
stopping inflation once it has begun [2]. Quantum physics can: accelerated 
expansion continually rips virtual infrared gravitons out of the vacuum [3] 
and these gravitons attract one another, thereby slowing inflation [3]. This 
is a very weak effect for GA <C 1, but a cumulative one, so inflation lasts a 
long time for no other reason than that gravity is a weak interaction [I] . 

This screening mechanism may be clear enough on the perturbative level 
but it has two frustrating features. The first is that, because inflationary par- 
ticle production is a 1-loop effect, the gravitational response to it is delayed 
until 2-loop order. The second frustration is that the 2-loop effect becomes 
unreliable just when it starts to get interesting. The effective coupling con- 
stant is GAHt and higher loops are insignificant as long as it is small. But 
all loops become comparable when GAHt becomes of order one, and the 
correct conclusion then is that perturbation theory breaks down. The break- 
down occurs not because any single graviton-graviton interaction gets strong 
but rather because there are so many of them. 

We believe it may be possible to derive a non-perturbative resummation 
technique by extending the stochastic method which Starobinsky devised for 

1 Here "large" means a A induced by a matter scale which can be as high as 10 18 GeV. 
Then, the value of the dimensionless coupling constant can be as high as GA ~ 10~ 4 
rather than the putative value of 10 -122 . 
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the same purpose in scalar potential models [H El C] • However, generalizing 
this technique to gravity is a difficult problem [HJ. This paper is part of 
an effort which is based on the idea of guessing the most cosmologically 
significant part of the effective field equations of quantum gravity. While 
there is no chance of guessing the full effective field equations, it might be 
possible to guess just enough to correctly describe the evolution of the scale 
factor a(t) for a homogeneous and isotropic geometry, using what we know 
from perturbation theory about how the back-reaction effect scales. Such 
simple cosmological models were recently constructed [9j[T0] and are reviewed 
in Section 2. Some basics of linearized gravitons needed in this work are the 
subject of Section 3. The possibility of getting an enhancement of high 
frequency gravity waves within this class of models was first investigated 
in [11]. In this paper we re- investigate this possibility using more accurate 
calculational techniques and we present our results in a way more appropriate 
for the needs of gravitational wave experiments. In Section 4 we review the 
enhancement mechanism while in Section 5 we derive an equation for the 
square of the magnitude of the mode functions and describe our improved 
evolution strategy. Our results are presented in Section 6. Their physical 
consequences and our concluding remarks comprise Sections 7 and 8. 



2 The Cosmological Model 

In a previous paper [9] we proposed a phenomenological model which can pro- 
vide evolution beyond perturbation theory. In one sentence, we constructed 
an effective conserved stress-energy tensor T^g] which modifies the gravita- 
tional equations of motion: E| 

Gfj, u = - -g^R = -Afiv + MGT^g] . (1) 

and which, we hope, contains the most cosmologically significant part of the 
full effective quantum gravitational equations. 

What form to guess for T^ u [g] was motivated by what we seek to do, and 
by what we know from perturbation theory. We seek to describe cosmology, 
which implies homogeneous and isotropic geometries. When specialized to 

2 Hellenic indices take on spacetime values while Latin indices take on space values. 
Our metric tensor g^ v has spacelike signature and our curvature tensor equals: R a g = 
r %, M + r % r % ~ (M v )- The initial Hubble parameter is 3H$ = A. 
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such a geometry the full effective stress tensor must take the perfect fluid 
form and we lose nothing by assuming that generally: 



T^[g] = (p + p)u fl u„ + pg^ . (2) 

The relation between p[g], p[g] and u^[g\ is heavily constrained by stress- 
energy conservation, but it is possible to specify one function for free. It 
turns out to be computationally simplest to take this free function to be 
the pressure [9j. We further require the pressure to be an ordinary function 
of some non-local scalar which grows like the number of e-foldings when 
specialized to de Sitter. If the pressure is to grow the way we know it does 
from perturbation theory [7], and to eventually end inflation, then a simple 
choice has the form [S]: 

p[g](x) = A 2 f[-GAX}(x) , X = i J2 , (3) 

where the function / grows without bound and satisfies: 

f[—GA X] = -GAX + 0[{GA) 2 ] , (4) 
and where the scalar d'Alembertian: 

□ = -/L^fVv^) , (5) 



-9 

is defined with retarded boundary conditions. The induced energy density 
p[g] and 4- velocity u^[g] are determined, up to their initial value data, from 
stress-energy conservation: 

= . (6) 

The 4-velocity was chosen to be timelike and normalized: 

gT u v ,u v = -\ u»u KU = . (7) 

The homogeneous and isotropic evolution [f| of this model - using a com- 
bination of numerical and analytical methods - revealed the following basic 



3 The line element in co- moving coordinates is ds 2 = —dt 2 + a 2 (t) dx-dx. In terms of the 
scale factor a, the Hubble parameter equals H(t) = da -1 and the deceleration parameter 
equals q(t) = -add -2 = -1 - H H~ 2 = -1 + e(t). 
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features: 

• After the onset and during the era of inflation, the source X(t) grows while 
the curvature scalar R(t) and Hubble parameter H(t) decrease. 

• Inflationary evolution dominates roughly until we reach a critical point 
X cr defined by: 

l-8nGAf[-GAX cr ] = . (8) 

• The epoch of inflation ends close to but before the universe evolves to 
the critical time. This is most directly seen from the deceleration parameter 
since initially q(t = 0) = — 1 while at criticality q{t = t cr ) = +|. 

• Oscillations in R(t) become significant as we approach the end of inflation; 
they are centered around R = 0, their frequency equals: 

u = GAH oy /72irft r , (9) 

where Hq is the constant inflationary Hubble parameter, and their envelope 
is linearly falling with time. 

• During the oscillations era, although there is net expansion, the oscillations 
of H(t) take it to small negative values for short time intervals - a feature 
conducive to rapid reheating; those of H(t) take it to positive values for 
about half the time; and, those of a(t) are centered around a linear increase 
with time. 

A novel feature of this class of models is the existence of an oscillatory 
regime of short duration which commences towards the very end of the infla- 
tionary era. During this period H(t) is positive about half the time, which 
represents a violation of the weak energy condition. Such a violation can- 
not occur in classically stable theories [12] but it can be driven by quantum 
effects of the type we seek to model without endangering stability [13] . 



3 Linearized Gravitons 

In terms of the full metric field gij{x) , the fluctuating graviton field h 
is defined as: 



TT, 



X 



9ij(t, 



a 



\t) 5 ij + V327iG /i^ T (t,x 



(10) 



4 In [9] , our analytical results were obtained for any function / satisfying (j4]) and growing 
without bound, our numerical results for the choice: f(x) — exp(x) — 1. 
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The free field expansion of the graviton field is: 

fff(t,x) = /^X:H^)^ k -%(k,A)a(k,A) + (c.c.)} , (11) 

where (c.c.) denotes complex conjugation, the polarizations ey(k, A) and op- 
erators a(k, A) obey: 

e ij (k ! A)e*.(k I A') = 5\\i , e«(k,A) = ^(k,A) = , (12) 
a(k,A),at( k ' 5 A')] = 5 XX , (2tt) 3 5 3 (k - k') , (13) 
and the mode functions u(t, k) satisfy: 



u(t, k) + 3H(t)u(t,k) + 



a 2 (t) 



it(t, fc) 







with the Wronskian associated with the two solutions of (tHJ) equaling: 

-3 



uu — uu 
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(14) 



(15) 



We shall be interested in the energy E(t, k) at time f of a mode with 
wavenumber k. The simplest way to derive this is to exploit the fact that the 
physical degrees of freedom of linearized gravitons have the same dynamics 
with those of a massless, minimally coupled scalar field (p(x). The scalar 
field Lagrangian density is: 



C(x) = - l -^~gcf v d^d uV = \a\t)<p 2 --Vtp-Vtp 



(16) 



The Langangian diagonalizes in momentum space: 

d 3 k 



L(t) 



d x C(x) 



(27T) 



|a 3 (t) k)' 2 



^a(t) A: 2 |^(t,k)' 2 



(17) 

so that any mode with wavenumber k evolves independently as a harmonic 
oscillator q(t) with time-dependent mass m(t) = a 3 (t) and angular frequency 
u(t) = ka- l (t): 



q{t) = u{t,k)A + u*{t,k)A ] , 
E S Ho{t) = ^a 3 (t)q 2 (t) + ^a(t)k 2 q(t) 



A, A f 



(18) 
(19) 



The analogous computation within the linearized graviton theory should only make 
an 0(1) change to the result. 
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At any instant t the minimum energy is E min (t, k) = |fca _1 (t). However 
since both the mass and angular frequency are time-dependent, the state 
with minimum energy at one time instant is not the state with minimum 
energy at another time instant; there is particle production as time evolves. 
The Bunch-Davies vacuum \Q) is the minimum energy state in the distant 
past and the expectation value of the energy operator (1191) in its presence 
equals: 

E(t, k) = ^a 3 (t) \ii(t,k)\ 2 + 1 k 2 a(t) \u(t,k)\ 2 . (20) 

A fair measure of the excess energy AE(t, k) acquired during time evolu- 
tion in any one wavenumber is obtained by subtracting the instantaneous 
minimum energy from f l20|) : 

AE(t,k) = (n\E(t,h)\Q) - E min (t,k) 

= ^a 3 (t)\u(t,k)\ 2 + h 2 a(t)\u(t,k)\ 2 - A . (21) 

4 The Enhancement Mechanism 

The oscillatory phase is a very distinctive feature of these models and in [TT] 
we investigated the possibility of gravitational wave enhancement due to its 
presence. There are two very plausible physical arguments that convinced us 
this is a worthwhile inquiry: 

- During the oscillations era the Hubble parameter H(t) changes sign and 
this, in turn, changes the sign of the "friction" term 2>Hu in the evolution 
equation (fT4l) obeyed by the mode functions u(t, k). For H(t) > this term 
tends to reduce \ii(t, k)\ whereas it tends to increase \u(t, k)\ when H(t) < 0. 
What happens to the magnitude \u(t, k)\ depends upon where u(t,k) is in 
its own oscillations when H(t) changes sign but the change from "friction" 
to "anti-friction" can clearly strengthen the amplitude in some cases. 

- The oscillations era is characterized by the frequency u given by (Q. Grav- 
itational waves of frequency close to u can resonate and their amplitude can 
increase. 

The first effort to evolve (TH|) through the oscillatory phase was done in 
[TT] . As expected, it is the near-horizon modes that experience enhancement: 
the natural time scale of their u(t, k) is close to the inverse of the oscillatory 
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frequency u and we get a significant resonance response, u When converted 
to current frequencies, the main conclusion of [TT] is the enhancement of 
gravitational waves with frequencies somewhat less than 10 10 Hz. In obtain- 
ing these results, however, certain assumptions were necessary since we do 
not possess exact forms for the two linearly independent solutions of f fl4|) 
during the oscillatory regime. Nor do we know which linear combination 
of these two solutions is the actual mode function as we do not know the 
linear combination coefficients. The latter are determined by knowledge of 
the initial conditions at criticality. Because the post-inflationary scale factor 
effectively describes an overall linear expansion on which the oscillations are 
superimposed [PJ, in [TT] we solved (1H|) for a linearly expanding a(t) - which 
can be done exactly - and then numerically superimposed the effect of the 
oscillations. We also had to make an "educated guess" regarding the initial 
conditions at criticality. 

In re-visiting the subject, we have developed a method - to be described 
in the next Section - which is considerably more accurate and, therefore, 
leads to robust conclusions. 



5 The Evolution Strategy 

• The Variable M(t, k) 

We wish to derive an equation for the quantity M(t, k): 

M(t,k) =u(t,k)u*(t,k) = \ u (t,k)\ 2 , (22) 
because it is directly related to the tensor power spectrum Aj[(t, k): 

Al(t, k) = 32nG ^ J a*x e"** (Q | h^(t, x) h^(t, 0) | fi> 

= l^.k 3 M(t,k) . (23) 

7T 

From the definition of M it follows that: 

M = iiu* + uu* , (24) 
M = uu* + 2uu* + uu* . (25) 



6 Here and throughout, super-, sub-, near- horizon is with respect to the modes fc cr whose 
first horizon crossing occured at criticality, when the transition from the inflationary to 
the oscillating era occured: fc cr = H(t CT ) a(t CT ). 

7 The actual mode function is the coefficient of the annihilation operator in the free 
field expansion of the graviton. 
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By using the fact that u satisfies (j!4p we conclude: 



2k 2 

M + 3HM H — M 

a 2 



By subtracting the square of (1151) from that of 
hand side of (1261) in terms of M and M: 



1 



uu 



4M 



M 2 + 



1 



2uii* . (26) 
we can express the right 

(27) 



and obtain the desired equation: 



M + 3HM + 



2k 2 



M 



2M 



M + 



(28) 



The goal is to find M(t, k) such that (T281) is obeyed. An exact solution is 
beyond our abilities but we can divide the full time evolution range into sepa- 
rate intervals and obtain reliable approximate expressions for M(t, k) within 
each of these. 

• The Evolution of M(t,k): Inflation 

- During the inflationary era, it makes sense to adopt a scheme that works 
accurately for any kind of mode and, at the same time, avoids numerical 
evolution for as long as possible. A method that seems optimal is the devel- 
opment of an asymptotic series expansion for M(t, k) in powers of H 2 a 2 k 2 : 



M{t,k) 



1 f, , , fHa\ 2 n . . /Ha 

_| 1 + aW (_) 



+ 



(29) 



Substituting the above in (1281) allows us to determine the leading coeffi- 
cients a(t), /3(t) of the series. The final form for the asymptotic expansion 
of M(t, k) becomes: 



M(t,k) 



2k a 2 



1 + (l 



+ 



+ 



7 3e 



21 n 3 o 



4 4 ) E 8H 2 



/Ha 



+ 



(30) 



As long as e does not get large, the series fl30l converges rapidly and we 
can use it to evolve all the way to within, say, 2 e-foldings before first horizon 



8 



crossing. § We shall, therefore, adopt this method and evolve very accurately: 
(i) any sub-horizon mode all the way to criticality, (ii) any near-horizon mode 
until, say, 2 e-foldings before criticality, and (iii) any super-horizon mode 
until, say, 2 e-foldings before first horizon crossing. Afterwards, in all cases 
equation (1281) is evolved numerically. 

The important cosmological parameters in the inflationary era are [9"j[TTj: 



a(t) 



-N 



H(t) ~ iu;W4iV+~ 



H(t) 
e(t) 



2H 2 



4iV + 



(31) 
(32) 

(33) 

(34) 



where N is the number of e-foldings before criticality. The initial conditions 
used in the numerical analysis are those inherited from (13TJj) at the appropriate 
time. 

• The Evolution of M(t, k) : Oscillations 

During the oscillatory era the important cosmological parameters are [9| ITT]: 



a(t) 
H(t) 

H(t) 
e(t) 



a CT C2 
to 



Ci + oj At + y/2 cos(u; At 
l-y/2 sm(u At + 4 



Ci + oj At 
~H 2 (t) - 



a/2 cos(w At + 



to- 



■\f2 cos(u; At + 



Ci + co At + y/2 cos(w At + 



y/2 cos(u At + (f)) 


Ci + u At + y/2 cos(u At + (j)) 




1 - y/2 sin(o; At + <p) 


2 



(35) 
(36) 

(37) 
(38) 



where At = t — t CT measures time with respect to criticality. In this regime 
we analyze equation (|28l) numerically. The parameters (0, C\, C%) in ( 1351 - 
1381) are chosen to match the outcome from the inflationary epoch ( 13TH341) at 



8 The error is in ignoring terms proportional to (-^r) and higher. Even when we reach 
2 e-foldings before first horizon crossing that is very small: = e~ 2 ) 6 = e~ 12 . 
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criticality, where N = and At = 0: 



= arcsm (^£!o) „ _I , (39 ) 

d = S in0- ^2 cos0 « 3 , (40) 

C 2 = \= « - . (41) 

Ci + v^cos^ 6 

• v4n Observable 

To connect with physical measurements, consider the excess energy AE(t, k) 
at time t of a mode with wavenumber k. It is given by equation (|2"T1) or, 
equivalent ly, by: 

a 3 M 2 r i 9 

AE{t,k) = + [2A;a 2 M-l] , (42) 

where we have used (127 1 122]) . We shall be interested in any excess energy AE 
acquired during time evolution through the oscillating regime. The resulting 
excess energy density Ap is: 

Perhaps of more relevance for gravity wave detectors is the amount of 
gravitational waves energy density Ap per wavenumber k, and divided by 
the critical density p cr : 

n gw (t,k) = — Ap(t,k) (44) 



dlnA; p cr rfln/c 

4G k 3 



AE(t,k) . (45) 



37r^ ow a3 (t) 

6 The Results 

We first define dimensionless variables: 

k a H 

t = u At , k = , ct ee — , ft ee — , .M = 2 A; a 2 M (46) 

LiJCLcr Q>rr U} 
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and re-express in terms of them the evolution equation 



d 2 M 



n 



dM 



+ H 2 (A- 2e) M + 



2k 



2 r 



dr 2 " dr 
the asymptotic series expansion (}2"9~]) : 



a- 1 



1 (dM 



2M \ dr 



M{t,k) = 1 4 
"9 21 



1 



e\ /Ha 



2' V 



ft" 



_4 8 4 
as well as the excess energy ( 142]) : 



7 3ex e 



4 4 8H 2 



Ha 



+ 



AS 



Co' 



aH 2 



and the observable f |45l ): 



M - 



i aM\2 



2H <9r 



+ 



ft" 



4a.M 



M - 1 



dlnfc gw 



„ i 4 1 /ft;\4 
Gw 2 x - — (-) 



X 



37T % 2 ^01' 

Ha 



AM 



A* 



1 dM 



2H dr 



M - 1 



We then discretize the time interval: 

t — y Tj = iAr , 



(47) 



(48) 



(49) 



(50) 



(51) 



and the resulting discretized evolution equation (|47|) determines Mi+2 in 
terms of Mi+i and A4,: 



M l+2 = 2M i+1 - Mi - Hi At (M i+1 - Mi) + H 2 At 2 (A - 2e) Mi 



(M i+ i - Mi) 2 2k 2 At 2 ( 1 

I i 



2Mi 



at 



M„ 



(52) 



The results of the combined evolution, using the asymptotic series ( 148]) until 
either 2 e-foldings before first horizon crossing or criticality - whichever 
comes first - and numerical integration of (152"]) thereafter, are presented in 



3 For super-horizon or near/sub-horizon modes, respectively. 
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Figures 1-14. It is important to note the following: 

- The conditions used for initializing the numerical integration at 2 e-foldings 
before first horizon crossing are provided by evaluating fj4*8|) at this point. 

- The conditions used for initializing the numerical integration at criticality 
are provided by matching the inflationary solution ( I3T1I341 with the oscilla- 
tory solution( l3l5ll3"8"j) at this point, so that the three parameters (<p, C\, C2) 
take on the values ()39H4ip . 

- The dimensionless wavenumber k ct which underwent first horizon crossing 
at t = t cr and is the wavenumber differentiating super-horizon from sub- 
horizon modes equals: 

2 

Her = U cr a cr = —j= « 0.38 . (53) 

In creating Figures 1-11 we have chosen 136 values of wavenumbers k ranging 
from 0.05 (super-horizon) to 10 (sub-horizon). 

- Inspection of Figure 14 makes evident the wild time dependence of the 
observable H 2 x (^-r i\w) • Note that had we not multiplied the observable 
by H 2 even wilder variations would occur when % passes through zero. Hence 
it is not clear to identify where the transition from the oscillatory to the 
radiation domination era took place. A reasonable assumption is that the 
transition occured when % > and e = 2. This determines the corresponding 
time t to equal: 

radiation =>- H > & e = 2 =>- r « 2n N osc . (54) 

- To study the effect of the duration of the oscillatory regime on the enhance- 
ment, we have displayed the results for N osc number of oscillation periods 
within the regime. The values of iV osc used range from A^ osc = - which 
corresponds to the usual transition from inflation to radiation domination - 
to Nosc = 10. Figures 1-11 present the results - for Ti > and e = 2 - in 
ascending order of N osc values. 

- The dimensionless factor Gu 2 is not included in the evaluation of (^-r ^gw) 
as seen in Figures 1-11, 14. 
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Figure 1 : Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, for zero periods of oscillations. 
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Figure 2: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginning of the first oscillation. 
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Figure 4: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginning of the third oscillation. 
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Figure 5: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginning of the fourth oscillation. 




Figure 6: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginnng of the fifth oscillation. 
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Figure 7: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginning of the sixth oscillation. 
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Figure 8: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginning of the seventh oscillation. 
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Figure 9: Fraction of the energy density per wavenumber divided by the critical density 
in gravity waves from our signal, at the beginning of the eighth oscillation. 
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Figure 10: Fraction of the energy density per wavenumber divided by the critical 
density in gravity waves from our signal, at beginning of ninth oscillation. 
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Figure 13: Time evolution of the Hubble parameter H during the oscillatory regime. 
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Figure 14: Time evolution of W 2 x ( . Q gw ) for k = 2 during the oscillatory regime. 
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7 Physical Consequences 



The results of the previous Section allow us to make the following remarks: 

• The existence of the enhancement effect is confirmed by our analysis. In 
Section 4 we argued, on physical grounds, that the effect is associated with 
sign changes of the Hubble parameter. This is explicitly seen in Figures 13- 
14 where there is a synchronization among the strongest enhancement peaks 
of the observable (Fig. 14) and sign changes of the Hubble parameter (Fig. 
13). Moreover, we see that the effect diminishes with time. 

• The far super-horizon modes left the horizon many e-foldings before criti- 
cality, their mode functions are "frozen" thereafter, and these modes are not 
affected much from the presence of the oscillating regime. 

• The near-horizon modes show the enhancement due to resonances close 
to the oscillatory era frequency u. Notice that as N osc increases the peak en- 
hancement magnitude increases and shifts towards higher values of k. Thus, 
it is the near / sub-horizon modes that feel the biggest enhancement. 

• The far sub-horizon modes show an ever- increasing "tail" with increasing 
k (Figs. 1-11). This enhancement is there and should be observed if very 
high frequency gravity waves become detectable in the future. It is present 
even when N osc = (Fig. 1) and, hence, it has nothing to do with the ex- 
istence or not of the oscillatory epoch. Furthermore, we can understand it 
without resorting to numerical results. When k ^> 1 the variable (1481) M. 
becomes essentially unity and the observable fl50l) simplifies considerably: 



Therefore, at any fixed time r, the value of the observable - being propor- 
tional to k 2 - will follow a parabola as k increases; this is explicitly seen in 
Figures 1-11. At any fixed k, the observable is proportional to a~ 2 and its 
value decreases accordingly as r increases; this is seen in Figure 14, albeit 
for k = 2. 

• The high-frequency "tail" discussed above will inevitably lead to ultravi- 
olet divergences. Ultimately it is the correct ultraviolet theory of quantum 
gravity that will have to address the issue. Nonetheless, an interesting ques- 
tion to answer - within the framework of ordinary perturbative quantum 
gravity - is which counterterms would absorb the ultraviolet divergences of 
our observable. 




M 



1 




(55) 
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As a first step, we use (HSjl to expand (loll in powers of (-^f ) 2 until we 
reach ultraviolet convergence. To make the connection with the available 
counterterms more direct, we also convert - using ()44j) - from the ratio Q 
to the excess energy density Ap. The result is: 



Su- 



rf 



d\n k 



Ap(t, k) 



UJ 



3 

2 £ 



3 9 e 



To make ( 1561) ultraviolet finite, two subtractions are needed: one to renor- 
malize the k 2 term and one to renormalize the constant term. As we shall 
see, the two counterterms which absorb the divergences are, respectively: 



Ad = 9l R 



-g 



AC, 



(57) 



These counterterms induce the following stress-energy tensor contributions: 

1 



AT 1 



AT 



20i 
2# 2 



2R 



]- gtlv R + 2(, 



R . 



(58) 
(59) 



We are interested in the energy density component of the stress-energy tensor 
for cosmologically relevant (FRW) spacetimes; in that case: 

Qgno 2 n 2 , 



FRW 
FRW 



AT 1 



AT; 



-144 #2 u 



- e - - e z + — 
2 4 2H 



(60) 
(61) 



and our assertion is established: ATq - given by ( 160]) - absorbs the quadrat- 
ically diverging order T-L 2 term in ( 156]) . while AT 2 - given by ( 15T]) - absorbs 
the order H 4 term in ( 156]) which diverges logarithmically, lj 

• The present value of the enhancement can be straightforwardly computed: 



d 



dlnk Qs 



^ ^matter ^ 



d\nk ngw . 



matter 



/ ^matter \ 



d 



0.3 x 10 



' \d In k 

' d 



n 



-3 



dlnk ngw . 



Guj 2 
4 x 10- 12 



(62) 
(63) 
(64) 



10 The value of g2 needed to subtract the order % 4 term in ([56]) agrees with that first 
found in 1974 by 't Hooft and Veltman [14]. 
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The value of the observable in the oscillatory regime, for given N osc and n, 
can be found in Figures 1-11. The remaining factor of about 10~ 15 in 
( I64p makes the enhancement effect very small and presently unobservable. 
The passage from fl63l) to (1611) is valid because after their entrance to the 
radiation era the gravitational waves behave like any other kind of radiation, 
and because during the radiation regime - unlike the matter regime - the 
product H 2 a 4 appearing in floTJ]) is constant. 

• The present frequency of the enhanced waves is given by: 

/now = ^— = ^ (— ) < l^Hz x e~-^ N < 10'Hz . (65) 

Z7T "now ^^T ^now 

For the estimate (|65|) we used the near-horizon value k ~ 1 as well as [11] : 

(_*E_) < e -63-iAiV ^ 10 -28 x e -lAN ^ (g6) 
"now 

u < 10 55 H now « 3.2 x 10 37 /fe , (67) 
where AN is the number of oscillatory e-foldings which we expect to be small. 



8 Epilogue 

From Figure dH it is evident that the signal is peaked at a narrow band of 
very high frequencies and is negligible at significantly different frequencies. It 
would be challenging to detect gravitational radiation at such high frequen- 
cies but detectors in that range have been proposed [IB] . As noted in the 
text, the phase of oscillations does not affect modes which experienced first 
horizon crossing more than a few e-foldings before the end of inflation. The 
wavelength of our effect is A = f^ 1 > 0.3m, whereas the smallest scale fea- 
ture which is currently observed in the cosmic microwave radiation is about 
10 22 m [IB]! Our model does not change either how matter couples to gravity 
or the propagation of linearized gravitons, so it has no effect on the spin- 
down rate of the binary pulsars. The gravity waves we predict will certainly 
distort how pulsar light propagates, but the short wavelength again seems 
to preclude a detectable effect. LIGO is not sensitive above frequencies of 
7000Hz, which is far too low. The situation is even worse with LISA's high 
frequency cutoff of 0.1Hz [IB] . 

11 As mentioned in Section 6, the displays of the observable ( ^ n fc O gw ) osc in all Figures 
therein lack an overall factor of Guj 2 . 
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